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ABSTRACT 

The Main Sequence (MS) of star-forming galaxies plays a fundamental role in driving galaxy evolution 
and in our efforts to understand it. However, different studies find significant differences in the 
normalization, slope and shape of the MS. These discrepancies arise mainly from the different selection 
criteria adopted to isolate star-forming galaxies, that may include or exclude galaxies with specific 
star formation rate (SFR) substantially below the MS value. To obviate this limitation of all current 
criteria, we propose an objective definition of the MS that does not rely at all on a pre-selection 
of star-forming galaxies. Constructing the 3D SFR-Mass-Number plot, the MS is then defined as 
the ridge line of the star-forming peak, as illustrated with various figures. The advantages of such 
definition are manifold. If generally adopted it will facilitate the inter-comparison of results from 
different groups using the same star formation rate (SFR) and stellar mass diagnostics, or to highlight 
the relative systematics of different diagnostics. All this could help understanding MS galaxies as 
systems in a quasi-steady state equilibrium and would also provide a more objective criterion for 
identifying quenching galaxies. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: high-redshift 


1. INTRODUCTION 

The stellar mass and SFR of galaxies are fundamen¬ 
tal quantities now being measured extensively, from low 
to the highest redshifts at which galaxies have been dis¬ 
covered. For star-forming (SF) galaxies the two quanti¬ 
ties are tightly corre l ated with each other, to the extent 
that, following Noeske et al.| (2007), such a correlation 
has got the designation of Main Sequence (MS) of star¬ 


forming ga l axies. In a series of seminal papers (Daddi 


eTaL]|2007 

Elbaz et al. 

2007 Noeske et ai.||2007) it was 

recognized that such tig 

nt correlation persists to at least 


redshift ~ 2 with nearly constant slope and dispersion 
compared to the c orrelation followed by SF galaxies in 
the local Universe (Brinchmann et aLj|2004). Many sub¬ 


sequent studies have followed, confirming the existence 


of a MS, all the way to at least z ~ 4 (Bernhard et 


al. 

2014; 

Karim et al. 2011; Kashino et al. 2013j |Mag- 

ne 

i et a 

201 

41 Pannella et al. 2009; Peng et al. 2010, 

hereafter P10 

Popesso et al. 2011 

2012 Rodighiero et 

al. 2010, 2011 

2014; Sargent et al. 

2012; Whitaker et al. 


sion and redshift evolution of the SFR— M* correlation 
can vary quite dramatically from one study to another, 
with the logarithmic slope of the relation ranging from 
0 .4 up to ~ 1 (see e.g. the compilation in Speagle et 


al. (2014)). In extreme cases, if galaxies are collected in 


TUlTmited fashion, no appreciable MS is recovered 
and the SFR runs flat with stellar mass (Erb et al.|2006 


Lee et al. 

2013 

Reddy et al. 

2006 

), a result dominatec 


with SFR above threshold are recovered (Rodighiero et 
aLpOljj ). 
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One reason why the MS of one group may differ from 
that of another stems primarily from how galaxies are se¬ 
lected for being star forming. One may select them by a 


(as in P10), or using the BzK two-color selection ( 

Daddi 

et al.12004 

20071 Pannella et al. 2009 

), or the rest-frame 

U V J selection (Whitaker et al. 12012 

2014 

Williams et 


al. 2009[) , or setting a minimum threshold for the spe¬ 


cific SFR = SFR /M* (Karim et al. 2011). Ultimately, 
all these criteria cut out galaxies with low sSFR, but 
the threshold differ from one criterion to another, and 
in a manner that may depend on mass, SFR, reddening 
and redshift. Criteria that tend to select as SF massive 
galaxies with low (but detectable) SFR will return a flat¬ 
ter MS compared to criteria that would not qualify as SF 
the same galaxies. Clearly, the use of different criteria 
makes less comparable to each other the corresponding 
results and may fuel sterile debates as to which of such 
criteria should be preferable. 

The importance of the MS comes from the fact that 
most stars in the Universe have formed in galaxies ly¬ 
ing around it within about a factor of ~2 in SFR. The 
mere existence of the MS, and its sharpness, indicate 
that there is order in nature, i.e., in the growth of galax¬ 
ies, as opposed to mere stochasticity. It also implies 
that both the mass and SFR of high redshift ga laxies 


must increase with time quas i-exponentially (Maraston 


et al. 2010 Renzini 2009; P10). The absolute value of the 


sSFR sets the clock of galaxy evolution (P10), setting the 
growth rate of individual galaxies and controlling their 
lifetime before they are quenched. The slope of the MS 
is directly connected to the (low-mass) slope of the mass 
function of star-forming galaxies, and if lower than unity 
would cause a runaway steepening of the mass fu nction 
if not contrasted by a tuned rate of merging (Peng et al. 


2014|. Moreover, much interest there is on the relation 
between the MS and the specific rate of mass increase of 
the dark matter halos hosting them, a proxy of the rate 


of gas inflow fueling the SF activity (Bouch det al.| 2010 
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Fig. 1. — The 3D SFR— M* relation for local galaxies in the SDSS database and 0.02 < z < 0.085. The third dimension is the number 
of galaxies in SFR-M bins. The drop towards lower masses is partly artificial, as no V/Vmax correction has been applied. This offers a 
clearer vision of the 3D structure, with the two prominent peaks, one for star-forming galaxies and one for the quenched ones. Notice the 
sharp ridge line of the SF peak, the extremely steep fall off in the number of galaxies, either way of the ridge line, the divide, which is then 
been taken as the definition of the Main Sequence of star-forming galaxies. On the North-West side of the divide one also encounter the 
starburst outliers, at 2 : ~ 2 include ~ 2% of galaxies and contribute ~ 10% of the global SFR (Rodighiero er al. 2011). The SE side of 
the divide is populated by a mixture of galaxies with lower SFR, with some being just in a temporary excursion below the MS band, while 
others are definitely on their way across the saddle, towards the peak of quenched galaxies. No V/Vmax correction was applied in order to 
have a better visibility of the two peaks. Data are from the SDSS database. 


Lilly et al. 2013 Peng & Maiolino 2014). Finally, and 
perhaps most importantly, the Si ll distance from the 
MS can be used to identify galaxies that do not qual¬ 
ify as belonging to the MS, such as starburst outliers on 
one side and, on the other side of the MS, those that 
have started the quenching process and, while still star 
forming, are in transition towards pure passive evolution. 

For all these reasons, a more objective definition of 
the MS of SF galaxies is in order: one which should not 
rely at all on such a pre-selection. In this paper we ac¬ 
cordingly propose a definition of the MS and suggest it 
should be generally adopted. Stellar masses and SFRs 
can be measured in many different ways, e.g., for the 
SFR using several different indicators through the whole 
electromagnetic spectrum, and the choice will depend on 
the available data and on the redshift of the sample. We 
emphasize that we are not proposing a universal way 
of measuring SFRs (and masses), our goal is instead to 
objectively define the MS once the measurements have 
been made , whatever measurement procedures were used. 
Still, such objective definition for the MS should also 
help quantify the relative systematics of different SFR 
and mass diagnostics. 

2. THE 3d SFR-MASS-NUMBER PLOTS 


of ~ 240, 000 galaxies for which SFRs have been esti¬ 
mated from their Ha flux and their stellar masses from 


We select galaxies from the Sloan Digital Sky Survey 
(SDSS) DR7 release ( |Abaza j i~an et al.||200 41 for lying at 
0.02 < 2 < 0.085 and having a reliable spectroscopic red- 


SED fits, following the same procedures as in Brinch- 


0.02 < 2 < 0.085 and having a reliable"spectroscopic : 
shift. Having excluded AGNs, this has provided a sample 


mann et al. (2004) and P10, with the exception that there 
SFRs were derived from SDSS DR4 release. Galaxies 
were weighted by the V/V mSiX method for masses below 
the completeness limit. This dataset is ideal for our pur¬ 
poses, because of its exquisite statistics, highly reliable 
redshifts and homogeneity in SFR and mass estimates. 
However, the same procedure can be repeated at higher 
redshifts as well, though other SFR diagnostics may have 
to be used. 

Figure 1 shows the 3D SFR-Mass-Number plot, con¬ 
sisting of the SFR— M* relation where the third dimen¬ 
sion gives the number of galaxies in fixed-size (0.2x0.2 
dex) SFR— M* bins. The two peaks correspond to ac¬ 
tively SF galaxies, on the West side, and quenched galax¬ 
ies, on the East side of the plot. Notice how sharp is the 
divide of the SF peak. Besides error effects, the saddle 
between the two peaks must include a mixture of galax¬ 
ies, with some being in a temporary excursion to a low 
SFR coming from the SF peak, others are in a tempo¬ 
rary excursion from the quenched peak, triggered by a 
minor gas-rich accretion event, whereas others are truly 
quenching galaxies on their final journey from the SF 
peak to the quenched galaxies depository. The number 
of quenched galaxies increases continuously from high to 
low redshifts, so across the valley the flow from the SF 
to the quenched peak must dominate. 
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In the 3D plot of Figure 2 the third dimension gives the 
product of the number of galaxies times their SFR, thus 
illustrating with extreme clarity where in the SFR— M* 
plane star formation is concentrated. Similarly, in Fig¬ 
ure 3 the third dimension now gives the product of the 
number of galaxies times their mass, hence illustrating 
where the stellar mass resides. 

In Figure 2 the quenched galaxies virtually disappear, 
as most of the star formation concentrates on the main 
peak. The modest bump noticeable South-East of the 
main SF peak is due to a low level of star formation 
still present in a fraction of quenched galaxies, although 
just upper limits to SFR might have been measured for 
many of them. For this reason, the shape of the quenched 
peak (especially on its low SFR side) should not be over¬ 
interpreted. A comparison of the two figures shows that 
the main peak has shifted from the star-forming to the 
quenched one, as indeed most stars in the local Universe 
are contained in quenched galaxies and bulges (Baldry et 


Figure 4 shows that the ridge line is straight linear up 
to the highest stellar masses in the sample, without a hint 


al. 


2004). 


A comparison of these three figures reveals that there 
is a fair number of galaxies in the valley between the two 
peaks (Figure 1) that indeed contribute to mass (Figure 
3) but not much to SFR (Figure 2). 


3. AN OBJECTIVE MAIN SEQUENCE DEFINITION 

As the importance of the MS stems from its dominant 
contribution to the global star formation, it is worth root¬ 
ing more deeply in this qualifying property its very def¬ 
inition. The sharpness of the SF peak in Figure 1 and 
Figure 2 offers indeed a quite natural definition of the 
MS as the ridge line of the SF peak, either in the 
3D number or in the 3D numberx SFR plots, as the two 
divides nearly overlap each other. These 3D figures of¬ 
fer unquestionable evidence for the existence of the MS 
as well as of its more natural, cogent definition. This 
ridge line coincides with the mode of the SFR distribu¬ 
tion at fixed M*. The median or the average SFR (at 
fixed M*) might be measured more accurately, but would 
depend on a pre-selection of star-forming galaxies. The 
ridge line, instead, does not require such pre-selection 3 . 


of flattening with increasing mass (see also Brinchmann 
|et al.|2004| A bending of the main sequence could be due 
to the growing fraction of the total mass being given by 
bulges already quen ched, hence cont ributing mass but no 
star formation (e.g., Whitaker et al.|2012|). We don’t see 
this effect in Figure 4 (up to ~ 10 fl IW Q ), even though 
bulges should be maximally developed at z ~ 0. Rather, 
the increasing bulge fraction with galaxy mass may be 
responsible for the global deviati on fr om ~ 1 of t he slope 
of the SFR— M* relation (Abramson et al.||2014). 

For galaxies away from the MS peak Figure 4 shows 
that quenched galaxies populate two distinct peaks, one 
at high mass and one a low mass, which result fr om the 
double-Schechter shape of their mass function (Baldry et| 
al. 2004} P10). The high-mass peak is well separated 


from the MS, with very few galaxies with intermediate 
SFRs. The low-mass peak is instead connected in a con¬ 
tinuous fashion to the MS, with many transition objects 
in between. This has a similar counterpart in the color- 
mass plot, where at low masses th e dist i nction betw een 
blue and red galaxies gets blurred (Taylor et al.12015}. In 


the frame of the P10 distinction between mass-quenched 
and environment quenched galaxies, the high-mass peak 
can be attributed to mass-quenching and the low-mass 
one to environment quenching. 

The sharp dichotomy between the high mass, quenched 
galaxies and those populating the MS arises from two ef¬ 
fects. The first, probably not the dominant one at this 
low redshift, is that the mass quenching process may be 
very fast, ensuring a rapid transition from the MS to the 
quenched galaxy repository. The second reason, that we 
believe is dominant at low redshifts, is that most massive 
galaxies were quenched a long time ago, so one expects 
very few to be in a transition phase now. Indeed, the 
stellar population properties of local massive early type 
galaxies (including ellipticals) indicate that in m ost cases 
q uenc hing took place some 10 Gyr ago (see e.g., Renzini 
2006, for a review). This interpretation is confirmed by 


the fact that virtually _all massive ETGs are already in 
place at z ~ 1 (Cimatti, Daddi & Renzini 2006), and 


to trace the MS (e.g., 

Karim et al. 

2011 

Pannella et al. 

sive galaxies to take place now. 

2009 

Rodighiero et al. 

2014) but it requires pre-selecting 

The situation is just the opposite at low-masses. The 


criterion. 

Figure 4 shows the projection of the the 3D surface in 
Figure 1 over the SFR— M* plane, along with contours 
for the number density of galaxies. Several notable fea¬ 
tures are emerging from this plot. The best straight-line 
fit to the ridge line is log(SFR) = (0.76 ± 0.01)log(M * 
/M 0 ) — 7.64 ± 0.02, with a slightly flatter slope than de¬ 
rived in P10 for u — g color selected blue galaxies in the 
same catalog, which was 0.9±0.01. This difference arises 
from the combination of two effects: the use of SDSS 
DR7 instead of the DR4 release; in addition, the u — g 
color cut included several galaxies with M* ~ 10 9 Mg 
and with SFR definitely below the ridge line, that we 
may interpret as quenching galaxies. 


environment quenching rate increases with decreasing 
redshift, following the growth of large scale structure 
overdensities, and therefore is maximum at redshift zero. 
Indeed, it is well know that the number density of 
quenched galaxies inc reases by a factor of ~ 2 since z ~ 1 


The ridge/ mode line was orig inally used by |Brinchmann et al.| 
rTl (see also |SaIim et al.|2007 1 to fit the SF sequence m SDSS, 


but all subsequent studies adopted a pre-selection of SF galaxies 
to define the MS. 


(e.g., Bell et al. 2004), which is largely due to the increase 
in the number of low-mass quenched galaxies, whereas 
the number of massive quenched galaxies is stable in this 
redshift range. 

To illustrate our proposal, this paper is limited to use 
data relative to the local Universe, but it is of great in¬ 
terests to follow the evolution of these 3D surfaces as a 
function of redshift, looking at how the twin peaks and 
the galaxy population of the ( “green”) valley between 
them are changing with time. Such surfaces can be used 
to trace the MS divide at the various redshifts, which 
should be feasible even if current galaxy samples at high 
redshifts are not as populous as those used here. This 
extension to higher redshifts is left to future works, but 


















































4 


Renzini & Peng 



Fig. 2.— Analog to Figure 1 where the third dimension gives the product of the number X SFR, hence showing where most of the star 
formation takes place. The ridge line of the star-forming peak is the new definition of the Main Sequence. The modest bump on the East 
side of the main peak is due to quenched galaxies, some indeed still supporting a low level of star formation while data can give only upper 
limits for many of them. The y max correction was applied for this plot. 



Fig. 3.— The same of Figure 2, where the third dimension now gives the product number X mass, hence showing where stellar mass 
is contained. The small SFR bump seen in Figure 2 has now exploded, as the majority of stars in the local Universe reside in quenched 
galaxies. The Umax correction was applied for this plot. 


here we can venture in making explicit a few expecta¬ 
tions. At redshift ~ 2 one expects to find a very dif¬ 
ferent pattern for quenched galaxies compared to what 
shown in Figure 4: mass quenching must have started at 


full steam (massive galaxies are growing very rapidly and 
very rapidly must be quenched). Instead, environment 
quenching has barely started, as overdensity contrasts 
are still small. One therefore expects the high-mass peak 




The Main Sequence of Star-Forming Galaxies 


5 



2.0x10 4 4.0x10 4 6.0x10 4 8.0x10 4 1.0x10 s 1.2x10 5 


Fig. 4.— The projection of the 3D surface shown in Figure 1 over the SFR—M* plane. Level contours for the number density of galaxies 
are shown, with colors ranging from blue to red as a function of number density. The V/V m &x correction was applied for this plot. 


of quenched galaxies to be already in place and grow¬ 
ing rapidly, whereas the low mass peak should be barely 
noticeable. The flow of galaxies across the green valley 
should be high at high masses, low at low masses. 

These considerations may help understanding why in 
XJV J- selected samples of SF galaxies | Whitaker et al.| 
(2012) find a flattening of the SFR— M* relation in high- 
ma ss gal axies at high redshifts, whereas | Whitaker et| 
ah] (2014) find a steepening of the relation for low-mass 


galaxies at low redshift. In the former case the flattening 
may be due to the inclusion of massive galaxies already 
in their mass-quenching phase, whereas the steepening 
at low redshifts is likely due to the inclusion of low-mass 
galaxies in their environment-quenching phase (such as 
those seen in Figure 4). 

Finally, Figure 5 shows a section of the twin peaks 
at log(M*) = 10.5, for both the number and the 
number x SFR distributions. The shape of the distribu¬ 
tion of the quenched peak is affected by the large num¬ 
ber of galaxies for which just SFR upper limits. More 
interesting is instead the distribution of the star-forming 
peak. The a of its best fit Gaussian is 0.3 dex, which 
comes from a com binati on of intrinsics spread and mea¬ 
surement errors (Salmi et al. 2012). We can notice two 
deviations from gaussianity in the wings of the distribu¬ 
tion. On the low SFR side the excess w.r.t. the best 
fit Gaussian is due to quenching galaxies, whereas at the 
opposite extreme the excess is likely due to starburst out¬ 
liers from the MS. Notice also that there is a ~ 0.2 dex 
shift in the peak of the pure number distribution and the 



log SFR m 0 /yr 


Fig. 5.— A cut through the twin peaks at log(M) = 10.5, 
in black showing the normalized probability distribution function 
(PDF) of the number of galaxies and in red the number x SFR 
distribution (solid lines). The red dashed line shows the best fit 
Gaussian distribution. 


number x SFR distribution, a shift that is independent of 
mass: the ridge lines of the two surfaces run parallel to 
each other with a 0.2 dex difference and one can choose 
one or the other as the MS divide. 




















6 


Renzini & Peng 


4. CONCLUSIONS 

We have proposed to objectively define the Main Se¬ 
quence of star-forming galaxies as the ridge line in the 
3D surface defined by the SFR-mass-number relation, or 
nearly equivalently in the surface in which the third di¬ 
mension is the product number xSFR. These surfaces can 
be constructed with very great statistical significance for 
the local sample of SDSS galaxies, and provide a vivid, 
cogent view of the reality of the MS as a major prop¬ 
erty of galaxy populations. Such definition can be ap¬ 
plied to samples galaxies at high redshifts as well, though 
with somewhat lesser statistical significance while dif¬ 
ferent SFR indicators may have to be used in different 
redshift ranges. However, with the advent of near-IR 
multiobject spectrographs SFRs from Ha can be derived 


up to z ~ 2.5, i.e., usi ng the same indicator as used here 


for local galaxies (e.g. 

Kashino et al. 2013 

Steidel et al. 

2014 

Wisnioski et al. 

2009). 


We emphasize that according to the new definition the 


MS of local, low—^ galaxies is indeed a straight line, 
with no sign of steepening at low masses or of flatten¬ 
ing at high masses, features that may emerge when pre¬ 
selecting star-forming galaxies before constructing the 
MS. The logarithmic slope of the SFR— M* relation is 
found to be 0.76 with the new definition of the MS, 
whereas it is 0.9 when using a u — g color pre-selection as 
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in P10, though part of the difference comes from using 
SDSS DR7 as opposed to DR4 data. 

A projection of the SFR— M * —Number relation over 
the SFR— M* plane reveals the existence of a number 
of low-mass galaxies with sub-MS SFRs, that we in¬ 
terpret as being undergoing environmental-quenching of 
their star formation, as expected in the phenomenolog¬ 
ical model of P10. At the opposite mass end, very 
few galaxies are now found in the course of their mass¬ 
quenching, as indeed the model predicts that most of 
mass-quenching should have taken place at higher red- 
shifts. 

In summary, we propose a definition for the MS which 
does not require a pre-selection of SF galaxies, which 
should facilitate the inter-comparison of results from dif¬ 
ferent groups, help understanding MS galaxies as systems 
in a quasi-steady state equilibrium and especially pro¬ 
vide a more objective criterion for identifying quenching 
galaxies. 
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